Particle Dispersion on Rapidly Folding Random Hetero-Polymers 
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We investigate the dynamics of a particle moving randomly along a disordered hetero-polymer subjected 
to rapid conformational changes which induce superdiffusive motion in chemical coordinates. We study the 
antagonistic interplay between the enhanced diffusion and the quenched disorder. The dispersion speed exhibits 
universal behavior independent of the folding statistics. On the other hand it is strongly affected by the structure 
of the disordered potential. The results may serve as a reference point for a number of translocation phenomena 
observed in biological cells, such as protein dynamics on DNA strands. 
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The study of random motion on complex structures is es- 
sential to the understanding of dispersion phenomena ob- 
served in numerous physical systems, ranging from epidemics 
spreading in complex networks and information transport in 
modern communication networks such as the internet 
In biological cells, the transport of macromolecules is ac- 
complished by a variety of translocation processes in which 
carrier molecules move along complex fibrous polymer net- 
works, e.g. myosin translocation on actin fibers |3] or trans- 
port on microtubules |4]. If the involved topologies are scale- 
free, diffusion is often anomalous, i.e. the mean square dis- 
placement of a particle violates the linear dependence on time 
(X 2 {t)) ~ P with < 7 ^ 1 i. Depending on the un- 
derlying microscopic dynamics, subdiffusive (7 < 1) as well 
as superdiffusive (7 > 1) behavior is observed. For instance, 
when a particle moves along a polymer in a complex fold- 
ing state, it may jump to a neighboring location in Euclidean 
space which is distant in chemical coordinates. Effectively, 
the particle moves superdiffusively along the chain 0,0] and 
performs a random walk known as a Levy flight. This mech- 
anism may explain fast target localization of regulatory pro- 
teins moving along DNA strands 1 8 ] . Levy flights have been 
observed in a variety of systems, ranging from chaotic sys- 
tems 01 and particle dispersion in turbulent flows full to for- 
aging animals | lH [l2il and climate changes (3- Levy flights 
are characterized by an exponent < fi < 2 which quanti- 
fies the degree of superdiffusion and is related to the heuristic 
dispersion relation X(t) ~ t 1 ^. When Levy flights success- 
fully mimic single trajectories, the associated stochastic evo- 
lution equations are no longer of the Fokker-Planck type but 
rather generalizations thereof which involve fractional differ- 
ential operators. Fractional models have contributed consid- 
erably to the understanding of these systems, in fact the terms 
fractional kinetics and fraction al dy namics have been coined 
to classify them l!4flaflra[l7lll8Lll9ll . Of particular interest 
are systems in which the cause for superdiffusive dispersion 
and the heterogeneity of the environment interact antagonisti- 
cally. 

In this Letter we introduce and investigate a model for su- 
perdiffusive particle dispersion on flexibly folding random 
hetero-polymers. We focus on the interplay between long 
range Levy type transitions due to folding and the quenched 




Figure 1: Random hopping along heterogeneous polymers. A: A 
particle (black disc) moving along the chemical axis x experiences 
a random potential V(x) associated with the random sequence of 
different types of monomers. When the chain is in a complex fold- 
ing state, locations that are distant along the chemical axis x may 
be close in Euclidean space (C,D). The folding topology is deter- 
mined by a connectivity matrix depicted in B where circular arcs 
indicate neighborhood in Euclidean space. C: A folding state with 
a characteristic mesoscopic scale a. D: A freely flexible chain with 
long range connections on all scales. Arrows indicate conformational 
change over time. 



random disorder caused by the heterogeneity of monomers of 
the chain. Based on simple assumptions on the hopping rate 
and configurational dynamics, we derive a fractional Fokker- 
Planck equation (FFPE) describing the motion of the particle 
along the polymer. We compute the relaxation properties as a 
function of the effective potential strength and the Levy expo- 
nent /i. We find that the dispersion speed depends consider- 
ably on fi, but becomes universal on larger spatial scales apart 
from a discontinuous change at fi = 2 (i.e. for ordinary dif- 
fusion). Furthermore, the relative concentration of monomers 
and thus the particular shape of the potential does not affect 
the ordinary diffusion process (fi — 2), but strongly affects all 
superdiffusive processes, a result we believe to be crucial for 
the understanding of transport phenomena in living cells. 

Consider the scenario depicted in Fig. [0 A particle is 
attached to a hetero-polymer and performs a random walk 
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along the chain. Let x denote the chemical coordinate with 
an inter-monomer spacing of a. The chain is flexible, and 
rapidly changing its conformational state defined by the Eu- 
clidean coordinate H(x) of each monomer. The heterogeneity 
of the chain is modeled by a potential V(x) which specifies 
the probability of the particle being attached to site x. In a 
thermally equilibrated system this probability is proportional 
to the Boltzmann factor exp [— j3V(x)]. The dynamics of the 
particle is governed by the rate w(x\y, t) of making a transi- 
tion y — ► x at time t. We assume that transitions occur only 
between monomer sites which are close in Euclidean space, 
i.e. when \R(y, t) — R(x, i) < a. We make the simplest pos- 
sible ansatz for this rate to take into account the requirements 
of Gibbs-Boltzmann statistics, the potential heterogeneity of 
the chain and the complexity of conformational states, 

w(x\y, t) = - e -^ v ^- v ^l 2 T{x, y; t). (1) 

T 

where the parameter r is the typical microscopic time constant 
and the function T(x, y; t) is defined by 



Y(x,y;i) 



1 if |R(a;,t) -R(y,i)| < 







otherwise. 



(2) 



T(x,y;t) reflects the dependence of transitions on the time 
dependent conformational state of the chain and is symmetric, 
i.e. T(x,y;t) — T(y,x;t). The function T(x, y; t) can be 
interpreted as a time dependent connectivity matrix (Fig.^J). 
The propagator p(x, t) of a particle initially (t = 0) at the 
origin evolves according to the master equation, 

dtp(x,t) = dy [w(x\y,t)p(y,t) - w(y\x,t) p(x,t)] (3) 



in which the rate is given by Eq. {0- The geometrical factor 
T(x, y; t) varies erratically and can be regarded as a stochastic 
process evolving on a time scale r g , which is generally differ- 
ent from the hopping time constant r. Averaging (denoted by 
[•]) over conformational states the dynamics reads d t [p] = 
[Cp] where the operator C is defined by the rhs of Eq. If 
conformational changes occur on smaller time scales than the 
hopping (r g <C t) we may substitute [Cp] ~ [C] [p], which 
represents a mean field approximation. In mean field, Eq. Q 
is given by 
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w(x\y,t)] = ± e -eW- v MV 3 pQx - y\), 



(4) 



where p(\x — y\) = [\R(x,t) — R(y,£)| < a] is the proba- 
bility that two given sites x and y are neighbors in Euclidean 
space. If the folding process is stationary, this probability is 
time independent, and due to translation invariance along the 
chain it is a decreasing function of distance in chemical space. 
The specific functional form of p(x) determines the asymp- 
totics of Eq. 10. Consider the situation depicted in Fig.QJ]., 
where the chain is knotted such that non-local transitions oc- 
cur on a typical scale a > a. On larger scales p(x) vanishes. 



In this case, a Kramers-Moyal expansion of the rhs of Eq. Q 
yields the FPE dtp — W'p + DAp, in which the diffusion 
coefficient is given by D ~ (a/a) 2 /t and the gradient force 
is determined by the potential V(x) along the chain. The sit- 
uation changes drastically for the type of chain sketched in 
Fig.^D. For a freely flexible chain the quantity p(x) follows 
an inverse power law with increasing chemical distance, i.e. 
p(x) ~ l/|x| 1+/ \ Typically p < 2 1 20] and thus p(x) lacks a 
well defined variance and consequently a typical scale in long 
range transitions. A particle moving along such a chain will 
behave superdiffusively and perform a Levy flight in chemi- 
cal coordinates. Inserting p(x) ~ l/|x| 1+A1 with < p < 2 
into Eq. l|4} and subsequently into Eq. the asymptotics is 
governed by a fractional Fokker-Planck equation (FFPE), 

D- X d t p = e~P y ' 2 A^ 2 e^' 2 p-p e^' 2 A^ 2 e"^/ 2 . 

(5) 

A detailed derivation is given in Ref. 11711 . Here, D is the 
generalized diffusion coefficient and the operator A^/ 2 is a 
generalization of the ordinary Laplacian, 



(A^/ 2 /) (x) = C, 



f(y) - m 



|1+M 



(6) 



with C M = 7r _1 r(l + p) sin(7r/i/2). In contrast to the ordi- 
nary Laplacian, A' 1 / 2 is a non-local, singular integral opera- 
tor, reflecting the superdiffusive behavior of the process. The 
boundary case p = 2 represents the limit of ordinary diffu- 
sion, i.e. Eq. reduces to an ordinary FPE. When the poten- 
tial vanishes, V = 0, Eq. (|5} becomes dtp = DA^' 2 p and 
is solved by the propagator of the symmetric Levy stable pro- 
cess of index p, i.e. p(x, t) = (Dt) 1 /! 1 (x/iDt) 1 ^) with 
L M (z) = (271-)- 1 J dk exp(ikz - \k\^) [17:]. In its general 
form Eq. (|5} describes the dynamics of Levy flights in exter- 
nal potentials obeying Gibbs-Boltzmann thermodynamics. 

In the following we investigate the relaxation properties 
of Eq. l|5} in random potentials V. Since the shape of the 
potential is determined by the ordering of different types of 
monomers along the chain, V(x) will be bounded and fluctu- 
ate about some average. Furthermore, it will generally possess 
a typical correlation length £. Without loss of generality we 
let (V) = and (V 2 ) = V^. The correlation length £ is de- 
fined by £ = V 2 Io°dyC(y) where C(y) = (V(y)V(O)) 
is the correlation function. The most straightforward way to 
incorporate these attributes into a model is by using Gaussian 
random phase potentials, 



V(x) 



1 

2^ 



dk<j>(k)< 



-i{kx+-8{k)) 



(7) 



which are defined by a set of random uncorrected phases 
$(fc) (with = — $(— k)) and the power spectrum S(k) = 
J dxe tkx C{x) (with cj)(k)4>* (k') = 2TiS{k)5{k - k')). The 
pdf associated with this choice of V(x) is Gaussian with zero 
mean and variance = (27T)" 1 / dfc S(k). Fig.0\- and|2^. 
show two realizations of random phase potentials, each one 
with a different power-spectrum (and correlation function). 
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Figure 2: A: A random phase potential with a power spectrum 
S(k) = 2£e(|fc - tt/2£|) where 6 is the Heaviside function. 
The correlation function C(x) = 2/tt s'm(nx/2£)/(x/£ : ) 2 decays 
in an oscillatory fashion (inset). The potential V(x) varies smoothly 
around zero. B: A potential with exponential power spectrum 
S(k) — 2£ exp [— 2£|&|/7r] and Lorentzian correlation function 
C(x) = [l + (7rs/2^) 2 ] . The potential V(x) shows more struc- 
ture on a finer scale. C: Copolymers with different relative concen- 
trations c of monomer types (gray and black). 




Figure 3: Ralaxation for various Levy exponents /i(inset) and ran- 
dom phase potentials. Panel A (B) corresponds to the potential in 
figure |2j\(B). lg[fc£] denotes the decadic logarithm. The gray lines 
indicate the limits G M (fc — * 0) for \i — 2 (upper line) and /i < 2 
(lower line), see Eq. 1 1 31 . Panel C depicts the generalized diffu- 
sion coefficient D M (e) for the potential in figure|2j\(B) and for three 
Levy exponents fi. The dashed and dotted lines are the results ob- 
tained from perturbation theory, i.e. (e) = 1 — 4e 2 (/i = 2) and 
D M (e) = 1 - e 2 (jU < 2). 



The relaxation properties are determined by the eigenvalue 
spectrum E(k) of the evolution operator C defined by the rhs 
of Eq. In order to compute the spectrum, the FFPE can 
be transformed by means of p(x, t) — e~^ vl - x ^ 2 ip(x,t). This 
yields a fractional Schrodinger equation with identical spec- 
tral properties, 



H 



—Tiip with 
D (-A^ /2 + U 



(8) 
(9) 



The operator H is symmetric and the effective potential U is 
related to the orig inal potential V by U(x) = e ev A^ 2 e~ ev , 
where v(x) = V/Vq is a rescaled potential of unit variance 
and e = /3Vo/2 = Vb/2/csT is the potential strength in units 
of k B T. 

For vanishing potential e = 0, we have U = and Ho = 
DA''' 2 which describes free superdiffusion when /i < 2. The 
spectrum of Hq is given by Eo(k) = Dk^ . The wave number 
k > defines the spatial scale of the corresponding mode. 
When a potential is present, the spectrum can be written as 
E(k) = Dfj,(k; e) fc M where D fl (k; e) quantifies the relaxation 
properties on scales ~ fc _1 with the unperturbed fc M -behavior 
as a reference. If D^{k\ e)/D < 1 the process relaxes more 
slowly compared to free superdiffusion. The spectrum E(k) 
can be obtained for weak potentials by perturbation theory. If 
e«l, the effective potential U can be treated as a small per- 
turbation, for U = 0(e). Up to second order in e the quantity 
D^(k\ e) reads 



where the effect on relaxation is provided by the function 

G fl (k) = ^fdqS(q)g fl (k/q) with (11) 



9n{z) 



(12) 



D^(k;e)/D 



(1 - 4e 2 G M (fc)) , 



(10) 



Figs. |5J\ and|3ji depict G^(k) as a function of k in units 
of the inverse correlation length £ _1 for the two types of ran- 
dom phase potentials defined in Fig.|2j\ and[2j}. The solid line 
depicts the limiting case of ordinary diffusion (p, = 2). The 
potential slows down the ordinary diffusion process (G M (k) > 
0) on scales larger than the correlation length, and speeds 
it up (G M (fc) < 0) on scales smaller than £. The function 
G M (fc) has a pronounced minimum at k w Moder- 
ately superdiffusive processes (fj, > 1) behave in a similar 
fashion, exhibiting the highest variation for k w On 
the other hand, G Al (fc) differs strongly for different [i in the 
asymptotic regime k <C Note also that on small scales 
(k > almost all processes relax faster than without the 
potential. In fact, G ft (k) < for 2 > // > fi c where 
/i c = 2 — In 3/ In 2 « 0.415. Surprisingly, this is no longer 
valid for strongly superdiffusive processes with /i < fi c . For 
instance, in the case fi = 0.2 (dashed line in Figs.|5J\ and[3ji) 
G^(k) is positive for k > implying that strongly su- 
perdiffusive processes are slowed down even on small scales. 
Comparing potential types, we see that the relaxation is differ- 
ent for each potential, but these differences become less im- 
portant in the asymptotic regime, which is governed by G M (k) 
as k — > 0. This limit can be computed from Eqs. il II I12L 
observing that L dqS(q) — tt and 172(2 — ► 0) = 1 and 
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Figure 4: The generalized diffusion coefficient D fl (e) for the 
copolymer potential, i.e. Eq. 1151 at three relative monomer concen- 
trations. The potentials (Fig.[2f ) either possess sparsely distributed 
peaks (c = 0.05), troughs (c = 0.95) or vary uniformly (c = 0.5). 
Dashed (/i = 2) and dotted (p — 1) lines represent perturbation the- 
oretic results. A: D M (e) for a Levy flight (n — 1) is different for each 
potential. B: -D M (e) for ordinary diffusion {p = 2) is independent of 
c, the curves coincide. 



9»<2{z -> 0) = 1/4, 



lim G M (fc) 



H < 2 
M = 2, 



(13) 



Hence, the asymptotic behavior is universal, with the excep- 
tion of the limiting case of ordinary diffusion, and is indepen- 
dent of properties of the potential. The range of validity of the 
limit dl3t . however, strongly depends on fi. The limit is not 
attained for marginal exponents (e.g. fi = 0.2 and 1.9) even 
on scales several orders of magnitude larger than the corre- 
lation length (Fig. |5J\ and |3j3). The above results are valid 
for small potential strengths e. For higher effective poten- 
tial strengths we investigate the asymptotics numerically. The 
quantity of interest is the normalized generalized diffusion co- 
efficient -D M (e) defined by 



D M (e) = lim D ^(k; e) / D . 



(14) 



In the perturbative regime Eqs. i ll Ob and H41 yield the univer- 
sal relation D M (e) = l-e 2 for^i < 2andL> AI (e) = l-4e 2 for 
ordinary diffusion. Figs.[3p and[3f) compare these results to 
those obtained numerically. Although the numerics deviates 
from perturbation theoretic predictions for greater potential 
strengths e, the universality still holds, i.e. the asymptotics 
(k — > 0) is independent of /i and of the statistical properties 
of the potential. The crucial property is the non-locality of the 
process (i.e. ft = 2 vs. fi ^ 2). Thus, as soon as the folding 
properties of the chain permit scale free transitions (p ^ 2), 
the behavior of D^(e) changes abruptly. 

The pdf of random phase potentials is symmetric with re- 
spect to the mean, i.e. a value V is as likely to occur as — V 
along the chain. For a number of hetero-polymers this as- 
sumption is inadequate. Consider the simple model copoly- 
mer depicted in Fig. |2jl!. The chain consists of a random ar- 
rangement of monomers, each one equipped with an intrinsic 
local potential parity w_ and v + , with — u_ = v+ > and an 
interaction range which we assume to be a Gaussian f{x—x n ) 
centered at the monomer site x n , 



with f{x) — exp[— x 2 /2a 2 } and a « a. The v n are randomly 
drawn from a pdf p(v) — c5(v — V-) + (1 — c)S(v — v + ). The 
relative concentration of low and high energy monomers is 
given by c and 1 — c, respectively. The parameter c determines 
the shape of the overall potential. When c<l/2(c>l/2) 
the potential consists of a series of localized peaks (troughs). 
Mean and variance of the potential are (V) = (1 — c) V- +c v+ 
and (CV - (V)) 2 ) = (1 - c)cSv 2 with 5v = v+ - v~. 
Figs.|4j\ and|4ji depict the results obtained for the generalized 
diffusion coefficient D^(e) on these types of copolymers for 
three values of c, each one representing one of the situations 
depicted in Fig.|2jr. The parameter Sv was chosen such that 
the variance is identical in all potentials. Although the value of 
-D p (e) in the weak potential regime (e <C 1) is consistent with 
the one observed in random phase potentials, for greater val- 
ues of e a striking deviation occurs. On one hand, the ordinary 
diffusion process {p = 2) is nearly insensitive to the shape of 
the potential, all functions Di (e) coincide. On the other hand, 
the superdiffusive process exhibits a more (less) pronounced 
decrease with increasing e when c = 0.95 (c = 0.05) as com- 
pared to an unbiased concentration of monomer types. 

The results reported in this Letter predict for superdiffusive 
behavior on folding polymers that the dispersion speed de- 
pends strongly on the specific arrangement of various types of 
monomers. This is in sharp contrast to the case of ordinary 
diffusion, which solely depends on magnitude variations of 
the potential. Therefore we expect that potential heterogeneity 
is an essential ingredient in superdiffusive translocation phe- 
nomena of proteins along biopolymers. 

D. B. thanks W. Noyes for interesting comments and dis- 
cussion. 



V[x) = 2_j v ™ f( x ~ X n)-> 



(15) 
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